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Abstract. A substantial fraction of the warps in spiral galaxies may result from bending instabilities if the disks are essentially 
self-gravitating. With N-body simulations, we show that galaxies with self-gravitating disks as thick as HI disks are subject to 
bending instabilities generating S-shaped, U-shaped or asymmetric warps. S-shaped warps persist during several rotations and 
keep the line of node straight. The warp amplitudes generated by bending instabilities remain however modest. Other factors 
must be invoked for extreme warped disks. However, bending instabilities can account for most of the cases reported in optical 
surveys, where the warp angle is generally less than 5°. This mode of warping is very sensitive to the disk flattening. It also 
constrains the fraction of dark matter distributed in the disk and in the dark halo. 
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1. Introduction 

Since the first observatio n of the Milky Way warp a t 
the end of the fifties JKerr & all Il956t iBurkel \l95Jk . 
half a century h as passed with out a clear explanation of 
this feature (see Binnev, 1992, for a review). The warp 
problematic comes essentially from its high frequency 
in spirals. In the local group, including t he Milk y Way, 
the three dominant spirals are warped (Roberts, 1966; 



Rogstad, Wri ght & Lockhari Il976t iNewton & E merson, 
19771) and among ne arby galaxies they are common (Sancisi, 
1976UBosmalll98ll). Statistics of warps in HI JBosmaHl99l( 
Richter & Sancisii Il994 iGarcia-Ruiz et all Il998h and in the 
optica l band ( Reshetnikov. Il998t iR eshetnikov & Combes. 
19991 Sanchez-Saavedra. Battan er & Floridol Il990t 

Sanchez-Saavedra et all 120031) reveal that more than half 



the spiral galaxies are warped and asymmetric. This implies 
that warps are either frequently or continuously generated, 
or persistent over several dynamical times. Disks are easily 
warped when subject to an external torque, genera ted for 
example by gravitational forces during interactions (Ouinn, 
1991; Herna uislll99 it [H uang & Carlberg .l 1997b or accretion 
events Jjiang & Binnevl 1 19991 iRevaz & Pfennig eJ l2001bt 



Lopez-Corredoira et al., 2002), or by magnetic forces in the 
gas ( Ba ttaner et all 11990). The challenge is to understand 
how isolated galaxies (in appearance) can conserve their warp 
with a relatively straight line of nodes (LON), avoiding the 
winding problem early noticed by Kahn & Woltier ( 1959). 
Solutions like normal modes of vibration (bending waves) 



have been proposed JLvnden-BellL 1 19651) . Unfortunately 
Hunter & Toomre ( 1969) showed that no discrete modes exist 
in disks with soft edges. If, in the presence of a dark halo, 
discrete modes can exist (Sparke, 1984; Sparke & Casertano, 
1988) they will be quickly damped by dynamical friction 



it- 
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Nelson & T remaine. fl99llDubinski & KuiikenLll995l) . 

However, as pointed out bv lSellwo od ( 1996), previous stud- 
ies have treated the case of a cold razor-thin disk. When tak- 
ing into account the velocity dispersions necessary to keep 
the disks stable against axisymmetric instability, disks may 
be subject to a bending instability (also named fire- hose in- 
stability). Bending instabilities were first studied by Toomre 
( 1966) in an infinite slab of finite thickness. The instability oc- 
curs when the centrifugal force of particles following a bend 
in the slab is higher than the gravitational restoring force of 
the slab itself. A formal description of t he instability can be 
found in iFridman & Polvachenkol i 1984 p.37). lArakil dl985h 
showed that an infinite slab is stable when the ratio of the 
vertical velocity dispersion cr 7 to the velocity d i spersio n in 
the plane <t n exceeds 0.2 93. ISellwood & Merrittl ( 1 1994 and 
Mer ritt & Sellwood (11994 have studied this instability in real- 
istic disks of counter-rotating stellar components. They found 
discrete bending modes m - and m — \ that c ould p ersist in 
disks even without counter-rotation. ISellw ood (1996) studied 
the growth of such modes in axisymm etric cases and found that 
they are strongly Landau damped ( IFridman & Polvach enko. 
1984). He also mentioned that some warm disks support long- 
lived axisymmetric flapping oscillations. 

Since the end of the sixties, most of the works on the warp 
problem have been built on the assumption of a disk embedded 
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in a massive non-rotating dark halo. However, there are now 
many observations that let us think that a substantial fraction of 
dark matter may be present in outer galactic dis ks, in the form 
ofrot ationally supported cold and clumpy gas ( Pfenniger et all 
fl994 : Pfenniger & Combes lTl994 . 

The presence of a weak but existing star formation rate 
far from the center of gala xies has been reported by differ- 
ent authors jN aeslund & Jo ersaeterl 1 19971 ISmith et all l2000t 
Cuilland re et all 12001). This star formation rate reveals that 
molecular gas must be present in abundance in the outer disk of 
galaxies. Mo reover, it has been pointed out long ago by Bosma 
dl978l Il98ll) that in samples of galaxies, the ratio between 
the dark matter and HI surface densit y is roughly c onstant 
well beyond the optical disk (see a lso ICarignan et all 1.1990; 
Broeils, 1992; Hoe kstra et all 120011) . These observations sug- 
gest that the dark matter follows the HI distribution along the 
disk. More recently, the study of the spiral pattern of NGC 2915 
has suggested that the very extended HI d isk in this galaxy is 
supported by a quasi self-gravitating disk (Bureau et al., 1999; 
Masset & Bureau, 2003). 

In the context of the warps, massive disks are interesting 
because the kinetic pressure support is then lower than the ro- 
tational support and su ch disks are in a regime where the Araki 
criterion (ArakiL I 19851) is relevant, thus bending instabilities are 
expected. 

In the first part, we examine the evolution of bending insta- 
bilities in different self-consistent galactic models, where the 
mass is dominated by a heavy disk. The models differ essen- 
tially in the thickness of the heavy disk. We show that between 
the very unstable thin disks and the stable thick disks, a regime 
exists where bending instabilities spontaneously generate per- 
sistent warps. In the second part, a semi-analytic analysis al- 
lows us to fix constrains on the validity of the previous results, 
when a dark halo is also taken into account. The values of the 
mass and flattening of the dark halo compatible with bending 
instabilities is determined. 



2. The galaxy model 

Our galactic mass model is composed of a bulge, an exponen- 
tial stellar disk and a heavy disk made of HI gas and dark matter 
proportional to it. 

1 . The bulge density profile is a flattened Plummer model: 



bulge 



R 2 z 2 Y 5/2 



where a - 1 kpc and b = 0.25 kpc. 
2. The exponential stellar disk takes the usual form: 



Pd(R,z) °c e 



■R/h s -wh. 



(1) 



(2) 



where the radial and vertical scale length are respectively 
H R =2.5 and//- = 0.25 kpc. 
3. To comply with the HI observations and slowly varying ro- 
tation curves, the heavy disk density profile is: 



Phd(R,Z)cc 



-Ru/R „-\z\/h z (R) 



R 



h z (R) 



<{r/r 




Fig. 1. Contribution of each component to the total rotation 
curve. The upper panel shows a classical representation of the 
rotation curve. In the bottom panel, the square of the circular 
velocity is plotted as a function of log/?. The latter representa- 
tion enhances the effective contribution of each component to 
the total rotation curve (Kalnais, 1999). 

For R > /?hd = 7 kpc, the heavy disk surface density de- 
creases as R l . It is smoothly truncated by the arccosine 
function at /?™ d ax = 35 kpc in order to keep an approximately 
constant rotation curve up to the disk edge. In agreement 
with observations showing t hat dark matter is not dominant 
at the center of galaxies JWeiner. Sellwood & Williams! 
2001; Bis santz et alll2003l) . the heavy disk surface density 
drops to zero at the center if R^ > 0. To be c onsistent 
with t he observed flaring of HI in the Milky Way (iBurtonl 
1992), the vertical scale height hAR) is radius dependent 
and is written as: 



h z (R) = h z0 e 



R/Rt 



(4) 



max\ 
hd j ■ 



(3) 



where h z o is the vertical scale height at the center and Rf the 
flaring radial scale-length. 

The respective masses of the three components are M\, - 
0.15 ■ 10" M for the bulge, M d = 0.46 • 10 11 M for the ex- 
ponential stellar disk, and Mhd - 1.61 • 10" M© for the heavy 
disk. With these values the rotation curve is approximately flat 
up to R - 35 kpc. Fig. ^shows the contribution of each com- 
ponent to the total rotation curve. The lack of HI or dark matter 
at the center of the galaxy implies a negative circular velocity 
squared contribution by the heavy disk component (see the bot- 
tom panel of Fig. ^. In the upper panel, where the norm of the 
circular velocity is traced as a function of radius, it produces a 
bump around 3 kpc. 

The initial vertical velocity dispersion cr z is found 
by satisfying the equilibrium solution of the stellar hy- 
drodynamic equation in cylindrical coordinates (see e.g., 
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Binn ev & Tremaina Il987l p. 199) separately for each of the 
bulge, stellar disk and heavy disk components, 



Pi °~ zi 



f 



dzpid z <&, 



(5) 



where the index i can label any component. These calculations 
are carried out on the polar coordinate potential solver of the 
Particle-Mesh TV-body code described in Pfenniger & Friedli 
( 1993). The ratio between the tangential velocity dispersion cr^ 
and the radial one cr R is derived from the epicycle approxima- 
tion (Binnev & Tremaine, 1987, p. 125): 



R- 3 d R (^ 3 5rO) z=o 



a\ 4Q 2 4R- 1 (d R ®) z=0 



(6) 



where k is the radial epicyclic frequency, and Q the rotation 
frequency, determined from the potential /^-derivatives as in- 
dicated in Eq. (|6j. Since in this approximation the radial os- 
cillations are decoupled from the vertical oscillations, the ra- 
tio crJcr R , which is precisely the Araki criterion, is not con- 
strained. This allows us to choose the radial velocity dispersion. 
We have used the following relation: 



cr 2 



o% 



= P Z 



■P 



R- 3 d R (R 3 d R ®\ =0 



(7) 



where v is the vertical epicycle frequency, determined from 
the potential z-derivatives as indicated in Eq. |7}. The fac- 
tor p is constant of order 1 and chosen to fix the Savronov- 
Toomre radial stability parameter Q - k ct r /3. 36 CL, where 
E = 2 f dzp(R, z) is the total surface density. 

The found velocity dispersions are then used to distribute 
the model particles in velocity space and to find the rotation 
velocity for each component with the Jeans equations. 

3. The simulations 

3.1. Parameters 

In order to explore the existence and evolution of bending in- 
stabilities, we have run a set of 9 models that differ in the thick- 
ness of the heavy disk, parameterized by Rj and hgo. Values are 
chosen such that the minimal ratio cr z /cr R in the disk ranges in 
an interval around Araki's limit . At equilibrium, cr\ is pro- 
portional to the thickness of the disk. From Eq. 0, the ratio 
<t z I<t r strongly depends on the second z-derivative of the po- 
tential. Thus, increasing the thickness of the disk increases the 
ratio cr z /cr R . Moreover, since cr z and cr R are linked by Eq. @ 
and since k and the surface density £ are almost independent 
of the vertical thickness, thinner disks are also more unstable 
with respect to the radial stability parameter Q. In Fig.|2] for 
each model, Q and cr z /o- R are traced as functions of the radius. 
With the parameter ft fixed to 0.77, all models are slightly sub- 
critical with respect to the generation of spiral density waves. 



1 Araki's criterion has been derived in the case of an infinite slab. 
However, it remains locally a good first order approximation when the 
rotation support exceeds the radial pressure support ( Sellwood, 1996). 




Fig. 2. The radial stability parameter Q (top) and the ratio 
<t z /<tr (bottom) as a function of the galactic radius. In both 
graphics, the curves correspond, from bottom to top, in the in- 
terval 10-20, to the models of increasing thickness, 01, 02, 07, 
03, 08, 04, 09, 05, 06, respectively. 
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Table 1. Parameters Rf and h z o for the different models. The 
respective flattening ( of the total potential is also indicated. 



Fig.|5]shows the distribution of the models in the plane <r z - 
cr z /o- R , at the radius R where <t z /<t r is minimal. Thin disks 
are found at lower left end while thick disks are found at the 
upper right end. The dotted line corresponds to Araki's limit. 
Parameters for each model, as well as (cr z /<7/{) min and the flat- 
tening ( of the potential are summarized in Tabled The poten- 
tial flattening is defined as the ratio of minor over major axis of 
the iso-potentials. 

The models are evolved with the Barnes-Hut treecode par- 
allelized to run on GRAVITOR, the 132-processor Beowulf 
cluster of the Geneva Observatory. To trust thin disk transverse 
instabilities it is important to not use an anisotropic Poisson 
solver such as a polar grid mesh, therefore, despite being much 
more CPU time consuming, a treecode approach is preferred. 
Each simulation contains 2 18 = 262 144 particles of equal 
mass. The softening length is 0.15kpc, and the time-step is 
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Fig. 4. Edge-on projections of the models 01, 
100x25kpc2 



02, 08, and 04. Times are indicated at the upper left. The box dimensions are 
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Fig. 3. Ratio <t z /<tr as a function of the vertical dispersion cr z 
atR- 15 kpc. The values are taken at the radius where ct : /o-r 
is minimum. The dotted line corresponds to Araki's limit. 



0.5 Myr. If not indicated, the length unit is the kpc, and the 
time unit is the Myr. 

3.2. Results 

As expected, the evolution of the models depends strongly on 
the thickness of the heavy disk. In Fig. [3] from left to right, 
models can be divided in four groups: 

1) Model 01 has a ratio <t z I<Tr of 0.18. It is very unstable. 
The bending instability occurs quickly and generates a transient 
asymmetric warp that extends up to z = 4 kpc at R — 35kpc 
(see Fig.0]a). The warp is associated with a spiral arm and its 
lifetime is about 1.5 Gyr. The ratio <t z /<tr is increased by the 
instability above Araki's limit, at about 0.4. After t - 3000, 
ct : /o-r is higher than 0.3 and the disk remains vertically stable. 
Face-on and edge-on projections of the model are displayed in 
Fig.0 

2) Models 02, 07 and 03 have still a ratio (t-Jctr well be- 
low the Araki limit. The bending instability occurs during the 
first 2 Gyr. An axisymmetric bowl mode (m = 0) grows dur- 
ing about 1 Gyr, before that cr z increases and stabilizes the 
disk. Fig. 0b) shows a spectacular U-shape warp of model 02 
at t — 2000. During the instability, a large bulge forms, pro- 
longed by a bar. The galaxies look like SBa types (see Fig. [5] 
t > 2200 Myr). In model 07, the mode is slightly less symmet- 
ric, due to the lopsidedness of the galaxy. A short-lived m = 1 



mode is present at the end of the instability, during the bar for- 
mation. 

3) The four models 08, 04, 09 and 05 develop S-shaped 
warped modes (m = 1). Except for model 05 which has a ra- 
tio <j z I<jr - 0.3 just above Araki's limit, all are unstable with 
respect to bending. In the case of model 08, the warp is long- 
lived and lasts more than 5.5 Gyr, corresponding to about 5 ro- 
tation times at R — 30 kpc (see Fig.[6j. In order to follow the 
evolution of the warped disk, we have divided it in a set of 18 
concentric rings from R - to 35 kpc. The inclination of each 
ring (warp angle) is determined by the direction of the major 
axis of the inertia tensor of the matter contained in the range 
[R - AR,R + AR], where AR - 1 kpc is the half width of the 
ring. The evolution of the war p is displayed in Fig.Qin form 
of Tip-LON diagrams (Bri22s, 1990). The diagrams represent 
the precession and nutation angle of each ring relative to the 
inner galactic disk. The outer dotted circle corresponds to an 
inclination of 3°. At t - 500, the disk is flat. The warp grows 
from t - 1000 to t — 3000 where it culminates (z - 2 kpc at 
R - 35 kpc). During the following 4 Gyr, the amplitude of the 
warp slowly decreases. At t — 6000 it is still clearly observ- 
able (see Fig.0]c)). Despite the differential rotation of the disk, 
the line of nodes (LON) remains quasi straight. However, in 
the o uter parts, it is slightly trailing, contradicting Briggs' third 
rule (Bri22s, 1990) 2 . The rotation period of the LON is about 
1 . 1 Gyr, which corresponds to the circular rotation period at 
the edge of the disk (R - 35 kpc). At about t - 6500 Myr, 
the warp disappears and gives place to a large bar. The insta- 
bility increases the vertical velocity dispersions by a factor 2 
at the maximum, around 10 kpc. At larger radius, the factor 
is 1.5 or less. Velocities in the plane are unaffected. At the 
top of Fig. [8] we compare the warp angle of the four models 
at the time where the warp is maximum. The warp amplitude 
of model 04 is the highest. At t — 3500 it reaches more than 
2 kpc above the plane defined by the central regions (see Fig.0] 
d)). However, its lifetime is slightly shorter than model 02 and 
the warp disappears after 3.5 Gyr at t — 5000 (see the bottom 
of Fig. |HJ. The slightly thicker model 09 and 05 also develop 
warps with similar amplitudes but their lifetimes are shorter 
than 2.5 Gyr. In all four models, the warp shape is character- 



The galactic rotation is counterclockwise. 
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Fig. 5. Evolution of model 01 and 02 between t - 200 and t - 4000, noted in the upper left corner. Columns 1 to 4 show the 
face-on and two edge-on (perpendicular) projections of model 01, Columns 5 to 8 correspond to model 02. The boxes dimensions 
are 50 x 50 kpc 2 
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Fig. 6. Evolution of model 08 and 04 between t - 300 and t - 6000. Columns 1 to 4 show the face-on and two edge-on 
(perpendicular) projections of model 08, Columns 5 to 8 correspond to model 04. The boxes dimensions are 50 x 50 kpc 2 
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Fig. 7. Evolution of the warp in model 08 from t - 1250 to 
t - 7000 Myr. The diagrams show the polar angles of the rings 
(tip-LON plots). The three dotted circles mark = 1,2,3°. 



ized by an inner fiat region that extends between 5 to lOkpc 

(Fig.Htop). 

4) Model 06 has a ratio cr z jcrR well above 0.3. In this 
model, no vertical instability occurs and the disk remains flat 
all along the simulation. 

4. The influence of a dark halo 

The previous simulations have been run in the absence of a 
dark halo, which has the advantage of allowing us to identify 
the warp origin uniquely to the self-gravity of the disk when 
Araki's instability criterion is met. We want now to check up to 
which point the disk self-gravity leads to bending instabilities 
when the disk is embedded in a conventional hot halo. 

Since the bending instability is closely related to Araki's 
criterion, we determine how the threshold is modified by dif- 
ferent flattenings and mass of a hot dark halo. We suppose a 
hot halo because these are less subject to their own instabili- 
ties. Thus, a rigid potential is a convenient approximation of 
the halo role, provided the velocity dispersions cr z , cr R in the 
disk are calculated consistently for the new local potential. 
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Fig. 8. Top: warp angle 6, as a function of radius, for the models 
08 (t = 3100), 04 (t = 3500), 09 (t = 2700), and 05 (t = 2800). 
Bottom: time evolution of the mean warp angle computed be- 
tween R = 25 and R - 35 kpc. 

4.1. Models with dark halos 

The models with dark halos are based on the six first models 
with Rf = 40 kpc. We simply add a quasi homogeneous thick 
mass distribution in the disk region, however bounded in ex- 
tension more abruptly than isothermal halos because we want 
to transfered mass from the heavy disk to the halo region of 
similar extent. Therefore we just adopt an inflated Plummer 
poten tial, the Miyamoto-Nagai potential (Mivamoto & Naaai, 
1975) with a radial scale-length H of 20 kpc: 

o- , GM - w 



^ 2 + [h - c + VFT?) 2 



where c is the vertical scale-height. Here the constants in the 
potential are expressed differently from the original Miyamoto- 
Nagai potential in order to better express the horizontal scale 
length with H. The halo flattening is parametrized by £/, = c/H. 
For R < H, (h is a very good approximation of the ratio of 
minor over major axis of the mass iso-densities. 

The mass fraction of dark matter that lies in the disk is 
parametrized by the factor /. We can thus write: 



Afhd = / Af dark , M h = (1 - f) M dark 

where Afdark - Afhd + Afh. The total mass is then: 



Af tr 



Af b + Af d + / Afdark + (1 - /) Afdark- 



(9) 



(10) 



4.2. Stability criterion 

In the previous section, we have shown that the Araki limit is 
a good stability criterion against bending oscillations if it is 
applied to the minimum ratio o- z I<t r along the galactic radius. 
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Since k 2 and v 2 are linear with respect to the mass, we can 
write (cf. Eq. Q): 



o~ 



-4=/r 



M b kI + M d k\ + M dark [/^ d + ( i_f )K 2 ] 



(ID 



M b v 2 + M d v 2 + M dark [/ v 2 d + (1 - f) vg] ' 

where the v, and /c, for each components are computed for a 
unit mass component. For the bulge, the exponential disk and 
the heavy disk, those values are obtained numerically with the 
Particle-Mesh potential solver. For the halo, these frequencies 
are calculated analytically. 

For each of the 6 models, the influence of the dark matter 
fraction distributed in the disk or in the halo as well as the halo 
flattening is explored by varying the parameters / and £/, in the 
range to 1. The results are plotted in Fig. [5] The Araki limit 
of 0.293 is traced with the solid bold line. The models that may 
be subject to bending instabilities (— < 0.293) are located at 
the bottom right of this line. On the contrary, stable models 
(^ > 0.293) are located at the upper left. 

In the absence of the halo (/ = 1), results corresponds to the 
previous Section. Models 01 to 04 are unstable and model 06 
is stable. Model 05 is slightly above the Araki limit. However, 
as already mentioned, it develops a weak instability. 

When the heavy disk is replaced by a spherical halo of 
equal mass (/ = 0), the potential is smoothed around the 
Z - plane and the vertical epicycle frequency v decreases. 
As k is weakly modified, the ratio — increases. Fig. [9] shows 
that the model is stable as long as <f/, is greater than 0.06. If 
£/, < 0.06, the halo plays the role of the disk, and allows the 
bending instability to develop. It is relevant to note that be- 
tween (h - 0.5 and 1, the halo flattening has little effect on the 
stability. If we consider dark halos with gj, > 0.5, which corre- 
sponds t o more realistic halos accord ing to cosmological sim- 
ulations (Dubinski & Carlbere, 1991) and weak lensing stud- 
ies (Hoekstra. Yee & Gladders, 2004), the maximal halo mass 
in the flatter disk case (h z o = 0.05) must not be heavier than 
0.8 M dark (= 0.58 M tot ). For thicker disks, only a halo with a 
mass comparable to the heavy disk or less (Mh < 0.36 M lot ) 
allows bending instabilities to develop. 

Some of our galaxy models can be used as first order rep- 
resentations of the Milky Way. If we assume that the Milky 
Way warp results from a bending instability and that the disk 
is most probably critical with respect to Araki's criterion, the 
solid lines in Fig. fallow us to fix an approximate constraint 
on the halo to disk mass ratio. With a vertical scale-height 
of about 320 pc at the solar radius (R Q - 8.5kpc), model 03 
has a similar thickness to the Milky Way HI (Burton, 1992). 
For any dark halo axis ratio larger than 0.3, the marginal sta- 
bility of the disk implies a dark halo to heavy disk mass ra- 
tio around 1 (Fig. |9j (c)). With this ratio, the surface density 
(\z\ < 700 pc) of each component in units of M /pc 2 are : 
E b = 0.4, Ea = 36.7, E hd = 48.0, E h = 3.7, for a total surface 
density of E to t = 88.8. The corresponding rotation velocity at 
the solar radius is 205 km s . 

5. Discussion 

In this work, the vertical stability of a galaxy model where a 
substantial fraction of the gravitational matter is contained in 



the disk has been explored, leading to several aspects discussed 
below. 



5.1. Bending instability vs. disk thickness 

The strength of the investigated scenario is to require only in- 
trinsic gravitational dynamics, therefore the reported effects 
are immediately produced whenever the stability threshold is 
reached, independently of spurious external torques produced 
by strong accretion events, galaxy interactions, a misaligned 
dark halo, or magnetic fields. In all the examined cases a suf- 
ficiently thin and massive disk produces a bending instability 
leading to a warp. 

Ou r res ults were already expected by iMerritt & Sellwood 
Jl994l) and ISellwood & Merrittl ill 994 in studies of counter- 
rotating disks, that m = and m - 1 bending modes sur- 
vive in disks without counter-rotation. Extremely thin disks 
with h Z Q < 150 pc (potential flattening f < 0.648) are strongly 
unstable and generate asymmetric transient warps, while thick 
disks with h z o > 550 pc (£ > 0.661, cr z /cr R > 0.3) are com- 
pletely stable. In between, m = (U-shaped warp) and m - 1 
(S-shaped warp) modes are excited, m - modes grow for 
about 1 Gyr and disappear when the ratio cr z /cr R exceeds 0.3. 
S-shaped warps appear for thicker disks {h z o - 250 - 450 pc, 
f = 0.652 - 0.658). They are characterized by a flat central 
disk and persistent straight LON, slightly trailing at large ra- 
dius. The LON rotates at the circular velocity period of the disk 
edge. The most interesting result is that in model 08, the warp 
is persistent during 5.5 Gyr, a substantial fraction of the galaxy 
life. It is not clear if this persistent mode may be akin to the 
flapping mode reported by Sellwood (1996) in axisymmetric 
disks. 

In our simulations, the two extreme disks (model 01 and 
model 06) differ in their density thickness by a factor of 1 1 . 
This translates to a potential flattening variation of only 2.5%. 
It is striking that observed galaxies have typical HI disk thick- 
ness in the range where bending instabilities can be triggered, 
because this transverse critical state would be similar to the 
disk radial state, which is also critical with respect to spiral 
arm formation. 

Indeed, in more realistic models taking into account the 
dissipational behavior of the gas, energy dissipation reduces 
the gas velocity dispersion, and gaseous disks can be expected 
to reach a regime marginally unstable with respect to bending 
instabilities. This is similar to the bar and spiral instabilities 
in the plane which heat radially a weakly dissipative disk to a 
marginal stability Q ~ 1 - 2 state. Thus, we must expect that in 
more realistic models including gas dissipation, bending insta- 
bilities will be repeatedly excited. 

Thus, a unifying picture of galactic disks emerges as dis- 
sipative systems maintained for a long time in a marginal sta- 
bility state, where the effects of energy dissipation are counter- 
balanced by dynamical instabilities in both the radial and verti- 
cal directions. Obviously star formation must also play an im- 
portant role in compensating energy dissipation. 
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(c) Model 03, /»zo = 0.25 



(d) Model 04, h z0 = 0.35 
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Fig. 9. Stability of the models with dark halos against the bending oscillations. The value of the minimum of — along the radius, 
is plotted as a function of the fraction of dark matter in the heavy disk / and of the halo flattening £/,. The solid line corresponds 
to the critical value of 0.293. The upper dashed line corresponds to the values of 0.4 while bottom dashed line corresponds to 0.2. 
Shaded regions correspond to unstable models. 
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5.2. Comparison with observed warps 

This idea of connecting gas dissipation to warps is sup- 
ported by observa tions. Warps a re more frequent among 
late type galaxies (Bo smal Il99ll) and they are completely 
absent among lenticula rs cha racterized by gas deficiency 
(Sanchez-Saavedraet al., 2003). 

It is also noteworthy that the present scenario predicts a pre- 
dominant occurrence of S-shaped warps, because in a slowly 
flattening disk due to energy dissipation, after the Araki insta- 
bility threshold is crossed the S-shaped warps occur before the 
U-shaped ones. 

As seen previously, only warps with modest amplitude 
are spontaneously generated, with warp angl e less than 5°. 
Extreme cases, like for example NGC4013 (Bottemaet al., 
1987), which has a 25° warp angle, seem to require stronger 
causes than an internal instability, such as tidal interactions. 
However, a s the w arp angle distribution in optical surveys 
peaks at 3° (Reshetnikov, 1998), this scenario may contribute 
to a large fraction of the S-shaped, the U-shaped, as well as the 
asymmetric warps observed in the optical. 

5.3. Relationship between warps and dark halos 

The semi-analytic analysis in Sect. |4]confirms our expectation 
that bending instabilities can occur as long as the self-gravity 
of the disk is locally dominant. If we accept that most galactic 
disks are thin and warped by bending instabilities, we obtain 
an additional constraint on the mass fraction and thickness of 
a dark spheroidal halo. For extreme thin disks, such as model 
02, the halo mass fraction with respect to the total mass (using 
Fig-Eland Eq. dl0» cannot exceed 0.47, for any halo flattening 
larger than ~ 0.3. For more typical disk thicknesses, such as in 
model 03, the maximum halo mass fraction drops to 0.33, for 
any halo flattening larger than ~ 0.3. The only way to obtain 
"classical" halos with a mass fraction of 0.9 is to decrease the 
halo flattening well below 0. 1 . 

It is interesting to note that if now we assume that the warps 
are produced by an almost marginal instability state, the above 
halo mass fraction limits become equalities for any flattening 
above ~ 0.3. Therefore the warps provide a sharp constraint on 
the halo relative mass, but a weak one on the halo flattening 
when the latter is above ~ 0.3. 

These results corroborate the recen t independent study of 
Tvurina, Khoperskov & Bizvaev (2004) where the halo mass of 
observed edge-on galaxies is estimated from the thickening of 
the disk due to the bending instabilities. 

6. Conclusion 

The most important result of this study is that long-lived galac- 
tic warps of observed amplitudes can be produced by internal 
bending instabilities, and that the dark matter to the total mass 
fraction within the warp radius and included in an extended 
hot halo of similar extent is limited to about a factor of 2 if 
warps result from a marginal bending instability. Recent studies 
linking the spiral arms and kinematic properties of disk galax- 
ies and the Milky-Way converge towards similar low values 



for dark halo masses within the HI disk range (Fuchs, 2003; 
Kalberla, 2003; Masset & Bureau! 1200 J) . 

Further simulations combining a heavy disk and a live dark 
halo including dissipative gas or star formation feedback will 
allow us to specify in more detail these first order dynamical 
effects. 
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